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ABSTRACT 



We use deep HST/ACS F555W and F814W photometry of resolved stars in 
the M81 Group dwarf irregular galaxy Ho II to study the hypothesis that the 
holes identified in the neutral ISM (HI) are created by stellar feedback. From 
the deep photometry, we construct color-magnitude diagrams (CMDs) and mea- 
sure the star formation histories (SFHs) for stars contained in HI holes from two 
independent holes catalogs, as well as select control fields, i.e., similar sized re- 
gions that span a range of HI column densities. The CMDs reveal young (< 200 
Myr) stellar populations inside all HI holes, which contain very few bright OB 
stars with ages less than 10 Myr, indicating they are not reliable tracers of HI 
hole locations while the recent SFHs confirm multiple episodes of star formation 
within most holes. Converting the recent SFHs into stellar feedback energies, we 
find that enough energy has been generated to have created all holes. However, 
the required energy is not always produced over a time scale that is less than 
the estimated kinematic age of the hole. A similar analysis of stars in the con- 
trol fields finds that the stellar populations of the control fields and HI holes are 
statistically indistinguishable. However, because we are only sensitive to holes 
~ 100 pc in diameter, we cannot tell if there are smaller holes inside the con- 
trol fields. The combination of the CMDs, recent SFHs, and locations of young 
stars shows that the stellar populations inside HI holes are not coherent, single- 
aged, stellar clusters, as previously suggested, but rather multi-age populations 
distributed across each hole. From a comparison of the modeled and observed in- 
tegrated magnitudes, and the locations and energetics of stars inside of HI holes, 
we propose a potential new model: a viable mechanism for creating the observed 
HI holes in Ho II is stellar feedback from multiple generations of SF spread out 
over tens or hundreds of Myr, and thus, the concept of an age for an HI hole is 
intrinsically ambiguous. For HI holes in the outer parts of Ho II, located beyond 
the HST/ACS coverage, we use Monte Carlo simulations of expected stellar pop- 
ulations to show that low level SF could provide the energy necessary to form 
these holes. Applying the same method to the SMC, we find that the holes that 
appear to be void of stars could have formed via stellar feedback from low level 
SF. We further find that Ha and 24/xm emission, tracers of the most recent star 
formation, do not correlate well with the positions of the HI holes. However, 
UV emission, which traces star formation over roughly the last 100 Myr, shows 
a much better correlation with the locations of the HI holes. 

Subject headings: galaxies: dwarf — galaxies: irregular — galaxies: individual 
(Ho II) — galaxies: ISM — stars: formation 
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Introduction 



The role of star formation (SF) in creating structure and shaping the interstellar medium 
(ISM) is an important but unresolved topic. Studies of the neutral hydrogen (HI) com- 
ponent in a wide variety of nearby galaxies, including our own, reveal numerous holes, 
shells, and bubbles (for simplicity, we will refer to al l types as h o les throughout the pa- 
per), which span a great range in size and age (e-K-i HeilesI 19791: Brinks &: Baiaia _ 19861: 
Puche et alJll992l:l0ev fc Massev"l995 l: iKim et aljll999l:IWalter fc Brinksl ll999l:lMuller et al. 



20031 : iRelafio et all l2007inChu 2008; Bagetakos et all 120091 ) . The origin of these ISM fea- 



tures has long been attribut ed to stellar feedba c k from massive stars, i.e., stellar winds and 
Tvpe II supernovae fSNe) flWeaver et al.l ll977l: ICash et al.l[l98ol : iMcCrav fc Kafatoslll987l : 
Tenorio-Tagle fc Bodenheimerlll988l . and references therein). Typical ages of these ISM fea- 
tures, m easured from the hole diameters and HI expansion velocities, range frona ~ 10^ — 10^ 
yr (e.g., Oev &: Massey 1995 : Walter fc Brinks 19991 : Hatzidimitriou et al. 2005 ). which sug- 
gests that even if massive O and B type stars, the primary contributors to stellar feedback, 
are no longer present inside the HI holes, we should still see remnant stellar populations 
(assuming a universal IMF). 

While the stellar feedback hypothesis of HI hole creation is appealing due to its intuitive 
simplicity, observational and theoretical studies do not provide much supporting evidence. 
The close proximity of dwarf galaxies in the Local Group (L G) has permitted detailed studies 
of correlation between stars and HI holes. In the LMC, iKim et al.l (119991 ) compared the 
positions of 103 giant or supergiant HI shells with Ha, a nd found a weaker correlation than 
expected given the young ages of the shells. Similarly, iBook et al.l (120081 ) find that while 
HcK traces sites of secondary triggered SF around superg iant HI shells in the LMC, i t is not 
exclusively associated with stars interior to the holes. iHatzidimitriou et al.l (12005 1) cross- 
correlated the locations of 509 HI holes in the SMC with known catalogs of OB associations, 
supergiants, Cepheids, Wolf-Rayet (WR) stars, SN remnants, and stars clusters, finding 59 
HI shells with no stellar component, and that holes without associated OB stars exceed 
those with them b y factor of ~ 1.5. Further, in the LG dwarf irregular galaxy (dl) IC 1613, 
Silich et al.l (120061 ) find that even the HI structures that do contain young stars do not have 
the observed properties consistent with the traditional stellar feedback theory. 

The studies of the origins of HI holes in galaxies beyond the LG show similarly mixed 
results when testing the stellar feedback hypothesis. In the pioneering study of M31, 



Brinks fc Bajajal (11986.) found a strong correlation between OB associations and HI shells 
smaller than ~ 300 pc. Correlating high resolution HI imaging of the M81 Group dl Holm- 
berg II with Ha, iPuche et al.l ( 119921 ) found st rong evidence for a stellar feedback origin to the 
51 HI holes. In a follow up study of Ho II, iRhode et al.l (119991 ) used integrated BVR pho- 
tometry to search for stars within the HI holes and found only 14% had a stellar component. 
A third study of Ho II concludes that Ha does not adequately trace the progenitor stellar 
populations; how ever, detected FUV emission associated with HI holes supp orts the stellar 



feedback model (IStewart et al.ll2000l ). In another M81 Group dl, IC 2574, iPasquali et al. 
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( 120081 ) use integrated colors from imaging from the Large Binocular Telescope to show that 
there is not a one to one correspondence betwee n stars and HI holes, while an evaluation 
of alternative models by ISanchez-Salcedd (120021 ) finds that HI holes in IC 2574 are likely 
formed by stellar feedback. 

Although a nurnber of studies of select individual holes have fo und associated young stars 



e.g.. lOey &: Masseylll995l : Ivan der Hulstlll996l : IWeisz et al.l 120091 ) . the overall evidence for a 



stellar feedback origin for all HI holes is not as compelling and a variety of alternate expla- 
nations have been proposed. For example, gamma ray bursts (GRBs) have sufficient energ y 
to create holes in the HI without leaving a remnant stellar cluster (e.g 
Simil arly, collisions between high velocity HI clouds and a galaxy (e.g., 
19861 ) or ram pressure stripping from the intergalactic medium (e.g.. 



Loeb fc Perna 1998). 



Tenorio-Tag le et al.l 



Bureau et al.. 2004 \ 



could provide the means to create observed ISM features. iRhode et al.l (|l999l ) suggest that 
modified energetics, either by an overestimation in the energies necessary to evacuate HI 
holes or by a non-standard IMF, could explain the observed discrepancies, although they 
acknowledge the latter is not likely. An alternate possibility is that grav itational instabil- 
ities and/or turbulen ce naturally give rise to structure i n the ISM (e.g.. lElmegreenI Il997l : 
Dib fc Burk"ertll200,5l ). The study of Constellation III bv iHarris fc Zaritskvl ( l2004l ) has cast 
doubt upon a number of these alternate theories and favors structure creation by numerous 
SF events. While it is likely that many of these mechanisms play a role in shaping the ISM 
at some level, quantifying the impact of stellar feedback on the ISM can give key insights 
into the feedback processes and self-regulation of SF. 

A further alternative is that the traditional theory linking stellar feedback to the cre- 
ation of HI holes is no t complete. Specifi cally, the most rigorous observational test of the 
stellar feedback model (IRhode et al.lll999l ) assumed that single age stellar clusters generated 
the SNe energy necessary to evacuate the HI holes. This is a reasonable assumption based 
on an extrapolation from the theore t ical models of the e xplosion and expansion of a single 
SN in the ISM (e.g., Chevaliei] 1974 : Weaver et al. 1977). In this scenario, SNe from a sin- 
gle generation of stars (i.e., a cluster) explode, releasing energy into to an isotropic ISM, 
thus creating the observed hole and quenching future SF, by sweeping up and heating of 
the gas interior to the hole. However, the study of spatially resolved SFHs in nearby dwarf 
galaxies has provided evid e nce that SF within ~ 500 Myr c an return to the s ame locations 
( iDohm-Palmer et al.l Il997l : iDohm-Palmer fc SkillmanI l2002l : IWeisz et al.l |2009| ) , demonstrat- 
ing the SF does not always suppress subsequent generations of SF. The creation of HI holes 
by multiple generations of SNe is a scenario that would likely have complicated dynamics, but 
preserves the intuitive simplicity linking the energy production of stars with the structures 
in the ISM. 

In this paper, we use HST/ACS imaging of Ho II to detect and study th e resol ved stellar 



populations within the HI holes indep endently catalogued bv iPuche et al. 



HI Nearby Galaxy Survey (THINGS; IWalter et all l2008l : iBagetakos et al 



dis, such as Ho II (roughly a solid body rotator in the inner 2 kpc; iPuche et al.lll992l ). are 



19921) and The 
2009li. Ne arbv 
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excellent environments for both the study of resolved stars and ISM features. Because the HI 
disks of dis are typically sohd body rotators, and thus lack sheer and differential rotation, HI 
holes and pote ntial progenitor stellar components t ravel in unison for time scales on order of 
~ 10® yr (e.g.. ISkillman et al.lll988l : ISkillmanlll996h . Further, the low metallicities of nearby 
dIs minimize internal extinction effects, which allow us to take a precise census of individual 
stars. 

From the resolved stellar populations of Ho II, we construct color-magnitude diagrams 
(CMDs) and measure star formation histories (SFHs) corresponding to catalogued HI holes 
that fall into our ACS field of view. Using the SFHs we can compute the energy associated 
with stellar feedback and compare it to the energy required to evacuate an HI hole, allowing 
us to quantify the efficiency of stellar feedback and compare it to predictions from models. 
We further compare the stellar populations of the HI holes to those of select control fields, 
similar in size to the HI holes, which span a wide range in HI surface density and do not 
overlap with catalogued HI holes. The resolved stellar populations a lso pe rmit us to directly 



test t he conclusions of previous studies of Ho II by iPuche et al.l (119921 ) and iRhode et al. 
( 119991 ). including comparisons to Ha, 24/im, and UV imaging of the same regions. Using 



the results of the stellar populations in the HI holes located within the ACS field of view, 
we use simulated CMDs to demonstrate the plausibility of HI hole creation due to stellar 
feedback in regions of low gas and stellar densities, such as the outer regions of galaxies. We 
then apply this method of analysis to the SMC and discuss how observations of apparently 
empty holes does not rule out the possibility of a stellar feedback origin. 



2. Observations and Photometry 



As part of a larger HST program aimed at imaging multiple dIs i n the M81 Grou p 
(GO-10605), we observed Ho II on 2006, December 30 using HST/ACS (iFord et al.l[l998h . 
To capture most of the optical galaxy, we used two ACS fields, observing each field for 
4660 sec in F555W (V) and 4660 sec in F814W (I), with a CR-split in order to reduce 
the impact of cosmic rays, and a dithering strategy to cover the chip gap. The images 
were processed with the HST pipeline and we performed photometry and artificial star tests 
using DOLPHOT, a version of HSTphot ( Dolphin 2000l ) optimized for ACS observation^. 
Photometry was carried out as part of a larger progr a m and the specific details, e.g., rejection 
criteria, completeness, can be found in lWeisz et al.l (120081 ). Following the photometric cuts, 
we merged the two fields into a single photometric list, appropriately removing duplicate 
stars, resulting in a total of 388,945 well measured stars for the combined fields. The deep 
photometry that is used in this paper has a limiting magnitude of My ~ 0, which allows us 



^Photometry and artificial star tests for t he HST/ACS i maging: of Ho II is available via the ACS Nearby 
Galaxy Survey Treasury program (ANGST: IDalcanton et al..,2009.) : http: / /www.nearbygalaxies.org_ 
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to see main sequence st ars (MS) up to ^ 300 Myr and Blue H elium Burning stars (BHeBs) 
up to ~ 1 Gyr in age (IGallart et al.ll2005l : IWeisz et al.ll2008l ). To quantify the accuracy of 
the photometry we ran artificial star tests on the entire galaxy and found 50% completeness 
limits of mi?555vi/ = 27.8 (Mp^ri ^w = +0.15) and mF8i4VF_ = 27.3 (Mpsuw = —0.15), assuming 



a distance modulus of 27.65 (iKarachentsev et al.l 120021 ). A number of HI holes in nearby 
galaxies have inferred kinematic ages (f rom HI observations) as old as a ~ 10* yr (e.g., 
Puche et al.lll992l : IWalter fc Brinkslll999D . If some of the older HI holes were created by SF 



on this time scale, only deep photometry of resolved stellar populations would probe faint 
enough limits to detect the remaining stars. 

The HI image we used in our analysis (panel (a) of Figure [1]) was observed with the 
VLA using the B, C, D configurat ions in 1990 and 1991 (the full detail s of the observation s 
are available in IPuche et"allll992l ). As part of the THINGS program flWalter et al.ll2008l ). 
the observations were reprocessed with their software pipeline, achieving a final spatial res- 
olution of ~ 6" and a velocity resolution of 2.6 km s~^. Both P92 and the THINGS team 
independently constructed catalogs of HI holes in Ho II by examining them in the position- 
velocity (pV) cuts and radius-velocity space in the HI data cube. Hole identification was 
done by eye using the following criteria as a guide: presence of the hole with a stationary 
center in multiple channel maps, sufficient contrast between the hole and its immediate sur- 
roundings, and the shape of the hole in ra dius- velocity space must be described by an ellipse 
flPuche et al.lll992l : iBagetakos et al.ll2009f ). Both catalogs contain the locations, diameters, 
expansion velocities of each of the holes (see Tables [1] — |4] and Figured]). Additionally, the 
THINGS catalog classifies each hole such that a type 1 hole has completely blown out of the 
disk of the galaxy, i.e, a break in the pV diagram, a type 2 hole only shows an approaching 
or receding side, and a type 3 hole has both sides present in pV space. While the P92 and 
THINGS holes catalogs are generally similar, the sizes, locations, and measured expansion 
velocities are not identical. For example, the THINGS catalog assigns an expansion velocity 
of 7 km s~^to any blown out hole, while P92 lists expansion velocities for each hole, and 
does not indicate a blow out. Because of the high stellar density, slight differences in the 
hole locations and size can lead to different stellar populations. We explore this further in 
U.3. 



With a distance of 3.40 Mpc, measured from the tip of the red giant branch (IKarachentsev et al 

0.05" 



20021 ). the high spatial resolution of HST/ACS (1 pixel = 0.05"; 1" ~ 16.5 pc) allows us to 
probe scales much smaller than the minimum detectable HI hole size (~ 100 pc). Treating 
each hole as circular and using the central coordinates and diameters from each hole catalog, 
we overlaid these apertures over the HST/ACS drizzled image and found that 23 of the 
51 holes from P92 and 19 of the 39 THINGS holes were within ACS coverage (Figure [2]). 
In the P92 catalog, a portion of holes 10, 21, 25, 34, and 48 are outside our ACS field of 
view. We included those that primarily overlapped with the ACS data (10, 21, and 48), but 
excluded the others from the sample (25 and 34). From the THINGS catalog, we excluded 
holes 18 and 38 for the same reason. For consistency we used the HI hole numbering systems 
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presented in the P92 and THINGS catalogs. Additionally, we selected nine control fields, 
regions in comparable size to the HI holes that do not overlap cataloged HI holes and span 
a range in HI column densities, and number them cl — c9 (Figures [Hand [2]). 



3. Stellar Components of HI holes in Ho II 



3.1. Luminous Main Sequence Stars Within the HI Holes 

From our photometry we constructed CMDs of the stars in each of the 23 holes in the 
P92 catalog, 19 holes in the THINGS catalog, and nine control fields (Figures [3] — [6] and 
Tables [H El & [5]) . Notably, we found that all holes in both catalogs contain a significant 
number of young stars; eno ugh to easily i denti fy CMD features by visual inspection and 
compare with the models of iMarigo et al.l (120081 ) (see Figure [7]). However, these CMDs are 
not consistent with the expected CMD for a single age stellar cluster. The holes also have 
red giant (RGB) and asymptotic giant branch (AGB) stars with ages < 1 Gyr, as expected 
given that t he older stars are f airly smoothly distributed throughout the optical galaxy (see 
Figure 11 in lWeisz et aDl2008h . 



Young stellar clusters 
mation of HI holes 



le.g.. 



e.g., OB associatiori s ) are traditionally associat ed with the for- 



Brinks fc Bajajal Il986l : iHatzidimitriou et al.l l2005l ) , and we tested 
this notion by making a census of the young est and brightest sta rs. We identified young 



stars based on the stellar evolution models of iMarigo et al.l (120081 ) . which suggest that the 
brightest 10 Myr old MS star at its turnoff has Mp^^^w ~ —5, or mp^^^w ~ 22.65 in Ho II. 
We found 302 MS stars with mp^^^w < 22.65 that also lay within holes in the P92 catalog. 
Further, restricting the color to be less than 0.05, i.e., mp^^^iY—mpgi^w < 0.05, (to exclude 
Red Helium Burning stars (RHeBs), bright foreground stars, and most BHeBs) we find 109 
stars that match our criteria. Of these 109 stars, 21 of them belong to hole 21, which is the 
largest and most populated hole with a diameter of ~ 2 kpc and ~ 40,000 stars in the CMD. 
Applying the above criteria for luminous MS stars to the remaining 21 holes, we found that 
all the holes have fewer than 10 bright MS stars and 20 holes have fewer than two bright MS 
stars. Similarly for the HI holes in the THINGS catalog, we found 347 stars with mp^^^w < 
22.65, and 104 stars with mp^^^w < 22.65 and m.F555W—'^F8uw < 0.05. Of these 104 stars, 
30 belong to the largest hole (number 17). The other 18 holes all have fewer than 10 young 
MS stars, and 13 have fewer than two young MS stars. For the stars located in the control 
fields (Figure E]), we found that 32 match the bright MS star criteria, with c7 containing 
21 of these stars, and six of the nine control fields containing two or less bright MS stars. 
The control field have more bright MS stars by a factor of ~ 1.5 when compared to the HI 
holes. We explore the impact of small number statistics, stochastic effects on the SFR, and 
production of massive stars in §5. 



From these basic calculations, it seems that the location of young OB type stars are not 
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a reliable tracer of the locations of HI holes, nor do they accurately represent the presence of 
young remnant stellar populations inside HI holes. We will discuss the spatial distribution 
of these stars further in 55.2. 



3.2. The SFHs of Stars Within the HI Holes 



Upper MS stars and BHeBs offer the opportunity to age date stellar populations more 
than 10 Myr old. To extract a more complete picture of the recent SF within the HI holes, we 
measured the SFHs for the stars inside each HI hole and compared the timing of SF activity to 
the inferred kinematic ages, i.e., the radius of the hole divided by its expansio n velocity. We 
derived the SFHs, from the CMDs, using SFH recovery code of [Dolphin fl2002h. This method 
constructs synthetic CMDs, using the stellar evolution models of iMarigo et al.l (120081 ). and 
compares them to the observed CMD using a maximum likelihood technique. To obtain this 
solution, we used a combination of fixed (e.g, binary fraction, IMF) and searchable (e.g., dis- 
tance, extinction) parameters. We allowed the program to search for the best fit metallicity 
per time bin, with the condition that the metallicity must monotonically increase with time 
toward the present. We found the mean metallicity in the most recent time bins from the 
SFHs to be consistent with the observed value of ~ 10% (IMiller fc Hodgdll996l ). which is 
best fit by Z=0.002 isochrones (Figure[7]). Throughout this paper, we use the following values 
when measuring SFHs or simulating CMDs: a standard power law IM F with x = —2 . 3 from 
0.1 to 100 Mq, a binary fraction = 0.35, the stellar evolution models o flMarigo et al.l ( 2008 ) 
a distance of 3.40 Mpc, a foreground extinction of Ap^^^w = 0.11 (ISchlegel el 
and 50% completeness limits of mp^^^w = 27.9 and mpsuw = 27.4 (IWeisz et al. 
quantify the errors in the SFHs, we added the systematic uncertainties from the isochrones 
and the statistical uncertainties from Monte Carlo tests in quadrature. Quantifying the best 



19981) 



2008h . To 



number of extensive studies (e.g.. 


Tosi et al. 


1989|; 


Tolstov & Saha 


1996|; 


Gallart et al. 


1996|; 


Miehell 


1997; 


Holtzman et al. 


1999 


Hernandez et al. 


199a 


DolDhii 


] 2OO2I: 


Ikuta & Arimoto 


2002 


: Yuk & Lee ,2007,). which explore the nuances of this method bevond the scope of this 



paper. Full details on the method we used in t his pa per to measure the SFHs and quantify 



the associated errors can be found in iDolphinI (120021 ). 



As an example we pr esent the observ ed, model, residual, and significance of the residual 
CMDs from the code of iDolphinl (l2002f ) for hole 23 in the THINGS catalog (Figure [8]). 
Visual inspection of the observed and model CMD show broad agreement of features such 
as the young MS, BHeBs, RHeBs, and the RGB. The difference between the model and 
the observed CMD can been seen in the lower two panels, where black points indicate more 
real stars than synthetic stars, and white points indicating more synthetic than real stars. 
The residual significance diagram (panel (d)), reveals no distinct patterns of black or white 
points, indicating that the model CMD describes the data quite well. 
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We chose to focus on SF over the past 200 Myr in each of the holes (Figures [9] — [TTl) 
and control fields (Figure [T2l) based on the inferred kinematic ages of the holes derived from 
both hole catalogs. The time resolution of our SFHs is ~ 10 Myr up to a look back time 
of 50 Myr and then ~ 25 Myr from 50 — 200 Myr. These choices provide an adequate look 
back time and time resolution to see potential SF events responsible for the creation of the 
HI holes. Note that kinematic ages for HI holes are traditionally considered an upper limit 
to th e age of the HI holes, as the expans ion velocities were most likely higher in the past 



e.g., iTenorio-Tagle fc Bodenheimeii 119881 ). We discuss the accuracy of inferred kinematic 



ages in more detail in §5. 

As a check for accuracy of SFHs between the two hole catalogs, we compared the SFHs 
and CMDs from overlapping holes. The CMD of P92 hole 44 significantly overlaps with 
THINGS holes 26 and 31. The CMD of P92 hole 44 has 7566 stars, where as the two 
THINGS holes have 4253 and 3190, respectively. In terms of the SFHs, the peak SFR for 
P92 hole 44 occurs in the 40 — 50 Myr bin with a peak values of 5.6x10"^ Mq yr~^. Both 
the THINGS holes how peak SFRs in the same time bin, with values of 3.0x10"^ Mq yr~^ 
and 3.2x10"^ Mq yr~^, or co-adding them, we find a peak SFR of ~ 6.2x10"^ Mq yr~^, 
consistent with the values for the single P92 hole. Similar tests were carried out for several 
other overlapping holes, and all were found to have consistent SFRs where expected. It 
is important to note that because the stellar density is not uniform in Ho II, even small 
differences in the location or radius of a hole, may result in significant differences in the 
stellar populations and SFHs. 

We searched for differences between the stellar populations of the control fields and the 
HI holes by comparing their cumulative SFHs (i.e., fraction of stars formed as a function 
of time) over the past 14.1 Gyr (lifetime) and 200 Myr (recent) by performing a KS test 
comparing the average cumulative SFHs for the holes and control fields (red lines in Figure 
[T3!) . As expected, the majority of the holes and control fields share a common old stellar 
population, i.e., at least 50% of the stars were formed greater than 6 Gyrs ago), and the 
fraction of stars formed in the last 6 Gyrs is also consistent between the hole and control 
field samples. Although the past 200 Myr is likely more influential when it comes to hole 
formation, again we see no discernible difference between the recent cumulative SFHs of the 
two samples. The KS test confirmed these impressions, indicating that, on both the lifetime 
and recent time scales, the stellar populations of the holes and control fields have the same 
parent distribution. From the perspective of the traditional stellar feedback hypothesis, this 
is an unanticipated result, as the theory would expect young stars associated with remnant 
clusters to be predominantly located inside the HI holes, and not control fields. This is 
discussed in more detail in 55. 
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Comparing the Energy from SF to the Energy Needed to Create an HI 

Hole 



Quantifying tlie energy generated from SF and the energy necessary to create an HI hole 
is of particular importance in understanding the role stellar feedback plays in creating HI 
holes. Fortunately, t he interaction of a single SN and the s urrounding ISM has been modeled 



in great detail (e.g., IChevalierlll974l : IWeaver et al.l 119771 : iHeilesI Il979l ). In addition, energy 



associated with SF is computed by the galaxy evolution modeling code STARBURST99 
( Leitherer et al.l 1999). STARBURST99 takes a SFH as input and uses stellar evolution 
models to produce predictions for spectrophotometric, chemical, and stellar energy output 
properties of galaxies that have had episodes of SF within the last 1 Gyr. For the partic- 
ular case of stellar feedback, STARBURST99 can compute the energy profiles (i.e., energy 
produced by stellar winds and SNe) using the SFHs for individual holes, which allows us 
to directly connect stellar feedback to the creation and evolution HI holes. In this section, 
we apply each of these methods to the HI holes in Ho II and their stellar populations, and 
compare them to expected stellar feedback efficiencies from a variety of hydrodynamical 
simulations. 



4.1. Energy Required to Evacuate an HI Hole 

The minir num amount of e nergy required to evacuate an HI hole is given by the single 
blast model of IChevalierl (119741 ): 



E 



Hole 



5.3 X 10^=^ A^.i2 ^1.4 



where EhoIc is in the units of 10^° ergs, no is the HI volume density (cm~^) at the mid- 
plane of the galaxy, d is the diameter of the HI hole in pc, and v is the expansion velocity in 
km s~^. While the expansion velocity and diameter of the hole are independently measured 
quantities, no is dependent upon the HI column density, Nhi, and scale height of the gas, 
h, in the following way: 



N, 



no 



HI 



2'Kh 



(2) 



The denominator, v27r/i is the effective thickness of the gas assuming a Gaussian dis- 
tribution, with h being the la scale height. 

One of the c hallen ges in determining Enoie is accurately measuring no. Both P92 and 
Bagetakos et al.l (120091 ) use a rotation curve to estimate the total mass in a gravitationally 
bound spherical system. This value was averaged over the volume of a sphere with a radius. 
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Rvmax ^iid multiplied by a f actor of two. Measuring ay = 6.8 km s ^(HI velocity dispersion), 
P92 found h = 625 pc, and iBagetakos eTaP (120091 ) find h = 600 pc. P92 used this method 
to compute the scale height of M31, and found it to be in agreement with the mass-model 
based method of iBrinks fc Burton! (119841 ). However, P92 acknowledged that wide-field, deep 
photometry of Ho II, would produce an accurate mass model of Ho II, and a more precise 
knowledge of h, hq, and Euoie- 

Fortunately, we had acce ss to high qual i ty 3.6/im imagin g of Ho II through the Local 
Volume Legacy Survey (LVL; iLee et al.ll2008l : iDale et al.ll2009l ) taken with the Spitzer Space 
Telescope, which allowed us to construct such a mass model. The LVL team provided a 
calibrated (including foreground star subtraction) 3.6/im image of Ho II. From this image, 
we used the ISOPHOT package in to perform surfa ce photometry, ac counting for 



inclination and position angle effects (z 
the dynamical center of the galaxy {a - 



= 47°, PA = 
08:19:18, 6 



177°; iPuche et al.lll992f ) starting at 



+70:42:37; J2000), which yielded a 
surface brightness profile a s a function of galactocentric radius (Figure [T^ . Following the 
method described in § 4 of lOh et al.l (120081 ) a nd assuming the 3. 6/im mass to light ratio, P* 
= 0.5 (comparable to other dwarf irregulars; IWalter et al.l 120081 ) . we converted the surface 
brightness profile into a surface mass density profile, E. 

Th e isothermal gas scale height equation (Equation 3) derived by KellmanI (jl970l ) (and 
used by lKim et al.lll999l for the LMC) relates the su rface mass density, S and t he HI velocity 
dispersion, ay, to the HI scale height, h. For Ho II. lBureau fc CarignanI (120021 ) derived ay as 
a function of galactrocentric radius, allowing us to compute the HI scale height as a function 
of radius. 



We now have two independent estimates of the gas scale height, namely the HI scale 
height of P92 and that inferred from surface photometry. We found the HI scale height 
derived by P92 to be larger by a factor of ~ 3 — 5 in the innermost 2 kpc of the galaxy, and 
in agreement at ~ 4 kpc. At radii larger than 2 kpc, h is in better agreement between the 
different methods. 

We computed new values of no, as a function of radius, using Equation (2) along with 
the newly calculated HI scale height and the values of Nhi measured by P92. Because the 
P92 Nhi values showed multiple large fluctuations in the inner region of the galaxy, likely due 
to the presence a number of HI holes, we interpolated over the regions of large fluctuations 
to ensure a smooth, continuous distribution. Effectively, the largest difference in the P92 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Asso- 
ciation of Universities for Research in Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 
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and new estimated values are due to the differences in the HI scales heights. The HI scale 
height derived from surface p hotometry effectively in creases the values of no by a factor of 



3 compared to the P92 and iBagetakos et al.l (120091 ) values, which in turn increases EhoI 



by approximately the same amount. 

It is important to note that the values of EhoIs are likely lower limits on the energy 
necessary to evacuate an HI hole. The values of no are the average HI volume density at 
a given radius. However, observationally, SF occurs at higher than average values of no, 
suggesting that to get a true sense of EhoIc we would need to measure a peak value of Nhi 
inside the HI holes. Although this is not observationally plausible, we can get a sense of 
how a larger values of Nhi would affect EhoIb- As an example, we examined control field 5, 
which is located near one of the highest regions of current SF in Ho II. From Equation (3) 
and the peak N^i measured inside c5, we find an increase by a factor of ~ 8 in the the peak 
value of no versus the average value. This translates into an increase in EhoIb by a factor of 
~ 10, suggesting that EhoIc values present here are likely lower limits. 



4.2. Energy Associated with Star Formation 

A basic criteria of the stellar feedback model is that the energy associated with SF must 
exceed the energy necessary to evacuate a given HI hole. To calculate the available energy 



due t o SF, we used the galaxy evolution modeling software STARBURST99 (iLeitherer et al. 



19991 ). We simulated a single instantaneous burst of SF with a fixed mass of 10 M0using 
STARBURST99 parameters that matched those described in §3.2, including a metallicity of 
Z = 0.002, a SN cutoff mass of 8 Mq, and a time sampling of 5 Myr. 

Because the energy output from STARBURST99 scales linearly with the input mass, we 
ran one simulation (with a stellar mass of 10^ M©), and then subsequently scaled it by the 
integrated stellar mass from the SFHs of stars inside each hole, sampled at 5 Myr intervals. 
The cumulative distribution of the energy generated from SF from the present (t = 0) to 200 
Myr for the stars inside each HI hole are shown in Figures [15] — |20l The inferred kinematic 
ages are denoted by the grey dashed line, and efficiencies (see §4.3) of 100%, 10%, and 1% 
are indicated by the red, green, and blue dashed lines, respectively. 

We have examined the efficiency over the inferred kinematic age for both catalogs as 
a function of galactro centric radius, hole radius, expansion velocity, and inferred kinematic 
age (Figure I^TI) . For both hole catalogs the clearest trend we found is that efficiency seems 
to increase with expansion velocity. However, this result is likely misleading, because of the 
differences in assigning expansion velocities. In the THINGS, catalog, if a hole was blown 
out, it was assigned an expansion velocity of 7 km s"^, whereas P92 attempted to measure 
the best fit expansion velocity for each hole. Thus, if we assigned each P92 hole with an 
expansion velocity lower than 7 km s~^, the value of 7 km s~^, there would be no trend 
in the P92 sample, consistent with the THINGS sample. Based on the THINGS catalog 
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classification (see §2) of liole types, we found that holes that have not blown out (type 2 or 
3), tend to have efficiencies higher than expected from models (see §4.3 for a description of 
the models), in some cases exceeding 100%. In contrast, holes that have blown out generally 
fall within the expected range. Again, this trend is likely due to the measurement techniques 
of the HI expansion velocities. In general, a blown out hole has stalled expansion, and its 
expansion velocity is not detectable above the random HI motions. Thus, the measurements 
of EhoIs are highly uncertain for blown out holes as are the efficiencies. The underlying 
implication is that the creation of an HI hole is a complex event that cannot be accurately 
represented by a simple kinematic age or efficiency, rather a more complex model is necessary 
(see §5 for more discussion). 



4.3. Comparing Energies and Feedback Efficiencies 



We then compared the inferred energy from SF over the kinematic age of each hole, 
EsFKA, to the energy needed to produce the observed properties of each hole, EhoIs- The 
ratio of EhoIb to Esfka allows us to calculate the a 'putative' stellar feedback efficiency, e, 
over the inferred kinematic age of each hole (see Tables [3] and Hj) . Models of SF including 
feedback and interactions with the ISM expect such efficiencies to range from ~ 1% — 20% 



(Theis et al. 


1992; 


Cole et al. 


1994; 



P92 holes and 6 of the 19 THINGS holes have putative efficiencies between 1 — 20%, in 
agreement with the models. Overall, a higher fraction of P92 holes are within the expected 
efficiency range than THINGS holes. While this may seem to be an unexpected result, 
this actually underlines the complexity of the connection between stars and HI holes. For 
example, P92 hole 42 and THINGS hole 27 are roughly the same size, but slightly different 
central coordinates, which results in the the THINGS hole having more stars in its CMD by 
a factor of ~ 2, in turn influencing the SFH and feedback efficiency. Of the holes that are 
not consistent with the expected range, one P92 hole (number 48) has an efficiency greater 
than 100%, while five THINGS holes (9, 12, 14, 23, 39) also have efficiencies greater than 
100% (Figure [2T1). These outliers suggest that perhaps energy from stellar feedback is not 
sufficient to have created these holes. 

For example, P92 hole 14 (Figure [T^ has an efficiency of ~ 1% where the inferred 
kinematic age intersects the cumulative energy profile. In comparison, P92 hole 30 (Figure 
[T6l) has a kinematic age efficiency of 39%, outside the theoretical range predicted by models. 
Further, ratios that require e > 100% indicate that there has not been enough energy for SF 
during the inferred kinematic age of each hole. However, Figures [TB] — [2D] show that there 
has been sufficient energy available during the most recent 200 Myr even in these extreme 
cases, suggesting that perhaps the inferred kinematic ages do not truly represent the ages 
of the holes, or that the concept of an HI hole age is intrinsically ambiguous for holes that 
contain multiple generations of stars. 
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To explore this possibility, instead of using the inferred kinematic age as a guide, we 
assumed a nominal 10% efficiency and recalculated the age of the hole to match the time 
scale needed to give sufficient energy input from SF. For example, while P92 hole 48 (Figure 
[T7I) has a kinematic age efficiency of 120%, the energy profile shows that over the past 200 
Myr, enough energy was generated from SF to have created the hole, if the hole is ~ 100 Myr 
old (assuming a nominal 10% efficiency). Similarly, for the other P92 and THINGS holes 
that have efficiencies outside the theoretical range, the same analysis shows that there is 
ample energy from SF to have created each of those holes in the past 200 Myr. We defined a 
'look back age' as the time at which the line of 10% efficiency intersects the energy profiles in 
Figures [15] — [201 We then compared the ratio of the look back age to the inferred kinematic 
age, and found that the inferred kinematic ages are generally lower than the look back times 
(Figure [221) . While perhaps not the most rigorous definition of the age of a hole, the look 
back age suggests that the inferred kinematic age may not accurately represent the true age 
of HI holes, even as an upper limit. We explore the kinematic age discrepancy in greater 
detail in §5. 



4.4. Connecting the Spatially Resolved Star Formation History to the HI 

Surface Density 

In addition to the measuring the SFR as a function of time, we can use spatial informa- 
tion from the resolved stellar populations to qualitatively compare the lo cation of SF episodes 
with HI column densities. Following a method of analysis similar to iDohm-Palmer et al. 



fll997h . we used BHeBs from the entire HST/ACS CMD of Ho II flWeisz et al.ll2008h and 
constructed movies of the spatially resolved recent SFH over the past ~ 200 Myr. The ini- 
tial masses of BHeBs range from ~ 2 — 20 Mq, giving them lifetimes between ~ 5 — 500 
Myr. More importantly, BHeBs have a one-to-one correspondence between luminosity and 
age. BHeBs of a single age therefore occupy a unique location on a CMD, making them 
superb tracers of recent SF. In contrast, multiple generations of MS stars can occupy the 
same location on the C MD, making age dating of most individual MS stars from the C MD 
virtually impossible (see iDohm-Palmer et al.lll997[ : iGallart et al.ll2005l ; IWeisz et al.ll2008[ . for 
a more in depth discussion of age dating stars from CMDs). 

By measuring its magnitude, we can assign a unique age to each BHeB, using the models 
of stellar evolution. We can then trace the location of SF as a function of age based on BHeB 
positions in the ACS image. While it is true that many you ng star clusters dissolve into the 
field on time scales of ~ 10 Myr (e.g., iLada fc Ladall2003l ). the stars can maintain spatial 
coherence for much longer. Empirically, CMD-based studies of the SMC and LMC find that 
star f orming structures ( such as associations ) can remain coherent for up to a couple hundred 
Myr flCieles et al.ll2008l : iBastian et al.ll2009l ). 



We were able to isolate the BHeBs on the CMD by visual inspection because Ho II has a 
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large number of BHeBs, and low foreground and internal reddening values. We constructed 
a smooth analytic function ma pping the magnitude s of the BHeB track to unique ages, using 
the stellar evolution models of Marigo et al. ( 20081 ). Applying this function to the observed 



BHeBs, we assigned each star an individual age, with an uncertainty of < 20%. We created 
density maps of the BHeBs by placing the stars into 5 Myr bins and converting their x and y 
coordinates to RA and Dec. We then smoothed the resulting images with a Gaussian kernel 
with a FWHM -15". 

To convert the BHeB density maps (stars pixel"^) to SFR/area (M© yr~^ kpc~^), we 
first divided each frame (i.e., 5 Myr binned density maps) by the sum of the pixel values to 
determine the fractional contribution of each pixel to the total stellar density in each frame. 
That is, if a normalized pixel had a value of it did not contribute to the SF in the frame, 
while a pixel with a value of 0.01 contributed to 1 % of the SF in th at frame. The density 



maps were then connected to the SFH of the field (IWeisz et al.ll2008l ). 



First, we re-binned the SFH with different time resolutions: 10 Myr for a look back time 
of ~ 100 Myr and 20 Myr for a look back time of ~ 200 Myr, in order to match the decreasing 
time resolution at older ages. We then interpolated each to a 5 Myr time sampling to match 
the temporal resolution of the stellar density maps. We then multiplied each stellar density 
map by the corresponding SFR and divided by the physical area of the ACS field of view, 
thus converting the images from units of stars pixel"^ to SFR/area. We then performed 
bilinear interpolation over time and space to ensure smooth and continuous movies. The 
result is two movies of the spatially resolved SFHs of Ho II with units of SFR/area, with 
look back times of ~ 100 and ~ 200 Myr, and a spatial resolution of ~ 15" (~ 250 pc). 
Figures [23] and [21] show select still frames of the two movies with HI contours overlaid. 

The general trends presented by the spatially resolved SFHs are that low level SF has 
dominated the galaxy in the most recent ~ 200 Myr, with only a substantial rise in the SF 
activity in the last ~ 40 Myr. The smoothness of the SFH implies that the HI holes were 
not produced by a series of intense episodes of SF, instead the holes were created due to low 
level SF with some variation above a constant value, which has steadily input energy into the 
ISM. As a comparison, we integrated the spatial resolved SFH of the region corresponding 
to THINGS hole 17 and found general agreement between the two. For example, both SFHs 
show a peak at ~ 40 Myr, and from the integration of the spatially resolved SFH we found 
the SFR at this time to be 3.0x10"^ Mq yr~^, whereas the SFH from Figure [S] shows a SFR 
of 4.0x10"^ Mq yr^^. While this difference of 30% may seem large, in fact, the systematic 
and statistical uncertainties from the two methods are ~ 20% from the temporal SFH and 
50% from the spatial SFH, which implies these are consistent values. 

These movies allow us to investigate the correlation between HI density and SF on 
different time scales. Within the last ~ 40 Myr, the SF is well correlated with high HI 
surface density, as has been seen in Ha, and UV imaging (see §5.2 for a discussion on Ha 
and UV imaging). On longer timescales, we note that the SF patterns do not adhere to the 
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locations of high HI column densities. In particular, there is a relatively high level of SF 
activity in the large HI void, corresponding to P92 hole 21 and THINGS hole 17, ~ 40 Myr 
ago, which has been preceded by low level activity for at least the past ~ 200 Myr. 



Discussion 



5.1. The Paradoxical Results of iRhode et al.l (119991 ) 



The study of Ho II by iRhode et al.l (Il999l ) is perhaps one of the most rigorous obser- 
vational tests of the stellar feedback model o f HI hole creation. Despite the appeal of the 
stellar feedback theory of HI hole formation, iRhode et al.l (Il999l ) found no compelling ev- 
idence supporting this hypothesis, and their conclusions seriously challenged the standard 
model of HI hole formation. The non-detection of stars in the majority of HI holes contra- 
dicted expec t ations for stellar feedback as a universal mechanism of hole creation, although 
Rhode et al.l (Il999l ) did not rule out stellar feedback as the source of small HI holes in the 
central regions of Ho II. 



Rhode et al.l (119991 ) searched for single age stellar clusters, but likely did not find them 
because the single age cluster model can overestimate the surface brightness of the remnant 
stars by a significant amount, for a variety of reasons we explore in the following sections. 
Further, using only integrated light, it can be difficult to distinguish between a single age 
stellar cluster and a mixed age stellar population, which can make it difficult to age date 
stars that are detected. As a com parative test, we computed the integrated magnitudes 
of each of the iRhode et al.l (119991 ) apert ures from th e HST /ACS CMDs. Our integrated 
photometry is in agreement with that of Rhode et al. ( 1999f). A histogram of the resulting 
integrated magnitudes reveals that the six regions where Rhode et al.l (1999) did detect stars 
are among the brightest in the sample (Figure [27]). The rest of the regions were classified as 
non-detections, despite containing significant numbers of young stars. 

Using the HST/ACS photometry of Ho II we construc ted CMDs of the apertures (panel 
(d) of Figure [1] and Table [6]) defined by lRhode et al.l (119991 ). Figures 125) and [26] unequivocally 
show well-populated CMDs of mixed age st ellar populat i ons, in cluding young MS stars and 
HeBs, which are likely the remnant stars [Rhode et al.l (Il999[ ) attempted to detect using 
integrated light. Our interpretation is that single age stellar po pulations do r iot pr oduce a 
majority of the HI holes in Ho II, confirming the conclusions of [Rhode et al.l (Il999l ). 



5.2. A New Model 



Instead of assuming that a single age cluster was responsible for the creation of HI 
holes, we explored the possibility that a mixed age stellar population producing SNe tens or 
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hundreds of Myr apart was responsible for the creation of HI holes. Specifically, we asked 
how a remnant mixed age stellar population would differ from a single age cluster in terms 
of the integrated light properties, location of stars inside HI holes, and energetics associated 
with this different mode of SF. 

To illustrate this, we generated synthetic CMDs, with the same parameters (e.g., IMF, 
binary fraction, distance, etc.) used to measure the SFHs in §2.2, and computed the mean 
expected integrated magnitudes of different stellar populations. To ensure that we were 
not in a regime where stochastic effects are important, we fixed the integrated mass of the 
simulations to be 10^ Mq, ~ 1% of the stellar mass of Ho II and at least a factor of 100 
greater than the stellar mass within any given hole. We limited the simulations to the most 
recent 100 Myr, the approximate age regime of the majority of HI holes in Ho II. The 
simulations were run 100 times, which produced an average statistical error of ~ 0.01 mags. 
Effects of completeness and photometric error are similarly ~ 0.01 mags. Neither source of 
error has a large impact on the results at magnitudes this bright. 

The first model assumed that all the stellar mass was formed in a ~ 5 Myr period 50 
Myr ago, with no subsequent SF. For this, 'burst' model, we found the mean integrated 
magnitude and color to be mp^^^w = 15.88 and mp^^^w ~ ^fsuw = 0.48. For comparison, 
we distributed the same integrated stellar mass with a constant SFR over the period of 50 
— 100 Myr and found mp^^^w = 16.07 and mp^^^w — mpguw = 0.51. 

We further considered a two component SFH model. We distributed 10^ Mq in stellar 
mass with a constant SFR from 50 — 90 Myr ago, no SF from 35 — 50 Myr ago, constant 
SFR from 25 — 35 Myr ago, and quiescence for the most recent 25 Myr. This multiple 
episode model produced stars with integrated properties of mp^^^w = 15.86 and mp^^^w — 
^F8uw = 0.54. Comparing the single age, constant, and multiple episode models of the 
same integrated stellar mass, we found they produced results that were indistinguishable in 
terms of integrated light. 

To illustrate this point empirically, we compared the stellar populations of P92 holes 
42 and 43, which both have integrated magnitudes of mp^^^w = 18.75, but differ in other 
properties. For example, hole 42 contains 1383 stars and has an inferred kinematic age of 17 
Myr, but hole 43 has 549 stars and an inferred kinematic age of 26 Myr. The SFHs of stars 
within the two holes differ significantly in the last 10 Myr. SF in hole 42 had been more 
prominent > 10 Myr ago, however within the last 10 Myr, hole 43 has a higher SFR, likely 
producing young bright stars that dominate the integrated light of the stellar population. 
These findings suggest that integrated light can only serve as an indicator of recent SF, but 
cannot give an accurate account of the SFH. 

The dissolution of stellar clusters into the field further comphcates integrated light 



observations and analysis of clusters. iLada fc Ladal (120031 ) note that the most common 



stellar clusters contain only dozens of stars and disperse into the field within ~ 10 Myr, 
which is the so-called 'infant mortality' scenario of cluster dissolution. Detailed studies 
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of structure of stellar groups in the SMC and LM C show that only ^ 3 — 7% of youri g 
stars survive in clusters past the age of ~ 10 Myr (IGieles et al.ll2008l : iBastian et al.ll2009l ). 
Applying this idea to Ho II, it suggests that the remnant stellar populations would not be 
highly clustered, but instead would be distributed across the HI hole, resulting in fainter 
integrated magnitudes than expected for clustered stars. To look for this effect we selected 
MS stars from the HST/ACS photometry by applying magnitude and color restrictions that 
isolated young MS stars (< 10 M yr (red) and < 75 M yr (blue)) and plotted their positions 
relative to the HI holes catalogs, iRhode et al.l (Il999l ) apertures, and control fields (Figure 
While there is some clustering of young stars, particularly those < 10 Myr in age, the 
clusters of stars are typically on the edges of the HI h oles where triggered secondary SF du e 
to the expansion of the hole is likely taking place (e.g.. iTenorio-Tagle &: Bodenheimerlll988l ). 
The majority of the MS stars < 75 Myr in age are not clustered, which makes the integrated 
light of these stellar populations appear fainter than had they been bound or appeared as 
unresolved stellar clusters. 



One of the most interesting cases presented by iRhode et al.l (119991 ) is that of P92 hole 
10. Based on its HI properties at least several dozen SNe would be necessary to have created 
hole 10. The remnant stars from a burst should be well above the detection criteria used by 



Rhode et al.l (119991 ) in a single cluster model, but none were detected. iRhode et al.l (119991 ) 
assumed a single remnant cluster would be in the center of the hole and used an aperture 
smaller than the hole to search for the stars. Comparing the CMD between the actual hole 
an d the aperture shows why the stars likely responsible for creating hole 10 were not detected 
by iRhode et al.l (119991 ). The CMD of the entire hole has ^ 4200 stars and an integrated 
magnitude of mp^^^w = 17.68, while the CMD of the Rhode et al.l (1999) hole 10 aperture 
contains 893 stars and has an integrated magnitude of mp^^^w = 19.75. Furthermore, the 
spatial distribution of MS stars < 75 Myr in age (Figure [28]) show that the photometric 
aperture of hole 10 has no stars, while the hole itself contains a number of young stars. It 
appears that the bright stars either did not form centrally in hole 10 or they have rapidly 
diffused out of the center, if they did form there. Although this hole was considered among 
the strongest evidence against the stellar feedback theory of HI hole creation, it now only 
seems to confirm that the single age cluster model of hole formation is not the appropriate 
mechanism for HI hole creation. 

If HI holes form via multiple generations of SNe, understanding the underlying physics 
becomes more complicated. We have to go beyond the simple single blast model of HI hole 
creation and consider the complex interplay between the local ISM and stellar feedback, 
which seems to be the do minant mechanism behind structure formation in the ISM (e.g., 
Harris fc Zaritskvl l2008h . iHeiled (Il990f ) considers the case of 'clusters of clusters' of SNe 
present within an HI hole, noting that sequential explosions would contribute to the me- 
chani cal energy exp a nding the hole as well as significantly alter the gas density inside the 
hole. lOey fc Clarkd (119971 ) also conclude that multiple OB associations would increase the 
mechanical power and time scale of energy input into expanding the shell. iRecchi fc Hensler 
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( 120061 ) present a model demonstrating how two separate SF events inside an HI hole would 
delay the onset of the hole refilling with cold gas with time scales from ~ 100 — 600 Myr 
depending on local conditions, indicating the multiple episodes of SF spread out over time 
can help extend the life of an HI hole. Projection effects could cause generations of stars 
along a line of site appear to be spatially co-incident within a given HI hole. However, this is 
likely a small effect because the disk-like nature of Ho II implies that the stellar scale height 
is much less than that of the HI, suggesting a reasonable assumption is that the young stars 
are primarily located in the mid-plane of the disk. 

Quantitative modeling of the impact of the complex scenarios presented by multiple 
generations of SF inside HI holes has yet to be carried out in detail. Presently, the single 
blast approximation and kinematic age estimations of HI holes are dependent upon quantities 
(e.g., expansion velocity, volume density) that likely have varied over the history of a given 
hole. Quantifying the precise impact of these parameters on HI hole formation and energetics 
is beyond the scope of this paper. For example, while why control fields and HI holes contain 
similar stellar populations, but vastly different HI column densities, is an open question. 
However, Ho II may not be the best galaxy to answer this question, because we are not 
sensitive to HI hole smaller than ~ 100 pc. Thus, if it is the case that larger holes are 
mergers of multiple smaller holes, we simply may not be sensitive to the small holes forming 
inside the control fields. Analysis similar to that presented in this paper using galaxies with 
smaller resolvable HI holes (e.g., LMC or SMC) along with computer modeling that accounts 
for the effects of an anisotropic ISM and multiple generations of SNe could lend insight into 
the fundamental physical processes that are responsible for shaping the ISM and driving HI 
hole formation. 



5.3. HI Holes in the Outer Regions of Ho II 



Of particular interest to iRhode et al.l (119991 ) are the presence of HI holes outside the 
the optical body of the galaxy. While they do not rule out the possibility that the inner 
holes could have formed from SF, the typical HI densities in the outer regions of Ho II are 
relatively low, which is not very conducive to SF (lKennicuttlll989l ). Howev er, the properties 
of ho l es in outer regions of galaxies do not differ frorn those in inner regions (IBrinks fc Bajaja 
19861 : IPuche et al.l Il992l : iHatzidimitriou et al.l l2005l : iBagetakos et al.l |2009|) , suggesting that 
the same underlying mechanism, namely stellar feedback, is responsible for HI holes in both 
environments. To this end, we used a Monte Carlo method to compute expected stellar 
populations in the outer holes of Ho II to see if their creation via stellar feedback is possible. 

From EhoIb we can calculate the SFR over the kinematic age, SFRka, for a typical 
outer HI hole and then use it to construct simulated CMD s to determine properti es of the 
expected stellar populations. We find that the 19 holes from IBagetakos et al.l (120091 ) located 
outside our ACS coverage have a mean diameter ~ 800 pc, an expansion velocity ~ 8 km s~^. 
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a kinematic age ~ 50 Myr, an HI volume density of 0.1 cm , and Enoie ~ 100 xlO erg. 
Assuming 10^^ erg per SN f|McCrav fc Kafatoslll987l) a nd that 10% of this energy goes into 



movi ng the ISM (iTheis et all Il992l : ICole et all I1994J : iPadoan et all Il997l : iThornton et al. 
19981 ). ~ 100 SN are needed to create an HI hole with these properties, equivalent to a SN 
rate over the kinematic age of ~ 10~^ SN yr~^. Assuming that the SNe come from stars 
with M > 8 Mq, we find SFRi^A ~ 2x10"^ Mq yr"^ for all stars, using our standard IMF 
(§3.2). This SFR implies an integrated stellar mass of ~ 10^ Mq over the kinematic age. 

With this SFR as input, we constructed simulated CMDs, via a Monte Carlo process 
using either a constant or a burst model. The input for the constant SFR model had a 
SFR = 2xl0~^ Mq yr~^ from 50 Myr to 100 Myr. This focuses on the brightest stars 
and most recent SF likely to have created the HI hole. The burst model assumed a burst 
of SF 50 Myr ago, producing 1000 Mq. We ran each simulation 100 times to account for 
the statistical fluctuations. For the constant model, we find that the average CMD of the 
remnant stars has ~ 7 stars and an integrated magnitude of mp^^ r^w = 24.00, assum ing a 



limiting photometric depth equal to the 50% completeness limit (jWeisz et al.ll2008l ). For 
the burst model, the average CMD of the remnant stars had ~ 8 stars and a integrated 
magnitude of mpH.^^w = 23. 7 5. Bo th cases results in integrated magnitudes fainter than the 



predictions by iRhode et al. fll999h as well as the hole and control fields presented in this 
paper. The remnant stellar populations in the outer holes are likely only detectable with 
deep photometry of resolved stars. 

Other common methods of searchin g for stellar populations inside HI holes include the 



use o f Ha, 24^m, and UV imaging (e.g.. |Puche et al.lll992l : iStewart et al.l 120001 : iMuUer et al. 
20031 ). UV traces a stellar population with initial masses > 3 Mq and time scales of ~ 100 
Myr, while Ha (and 24/im in more dusty environments) t races ionizing p hotons from MS 
stars with initial masses of > 20 Mq and ages < 10 Myr (see lKennicuttlll998l . for a full review 
of different SFR indicators and their properties). A number of studies aimed at testing the 
stellar feedback theory of hole creation via this type of imaging have have done so with a 
mixed degree of success flPuche et allll992l;lKim et al.lll999l: iRhode et al.lll999l: IStewart et al.l 



2000l : iMuUer et al1l2003l : ICannon et aDl2005l : iBook et al.l[2008h . Given the low SFRs we have 
derived (10""^ - 10"^ Mq yr"^) within many of the HI holes, it is not entirely surprising that 
these imaging techniques have yielded such mixed results. Simulations have shown that for 
SFRs below ~ 10~^ M^n yr~^ stochastic p rocesses become iinportant and Ha is no longer a 



reliable tracer of SF fiTremonti et al.ll2007l : iThilker et al.ll2007l : iMeurer et al.ll2009f ) . Ahhough 



a similar study in the UV has yet to be conducted, given the small number of stars expected 



i n any given CMD, it is likely that the UV has a similar cut-off SFR. Indeed, IStewart et al. 



( I2OOOI ) used FUV imaging of Ho II to show that there is some agreement between HI hole 
locations and bright UV regions, but did not find a good one to one correlation. 

We can test for UV-bright stellar populations within the out er holes using data from the 
Local Volume Legacy sample ( Lee et al.ll2008l : iDale et al.ll2009l ). w hich contains FUV an d 



NUV imaging originating from the GALEX Nearby Galaxies Survey (iGil de Paz et al.ll2007l ) 
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and the Hct and 24/xm images from the Spitzer Inf rared Nearby Galax ies Survey (SINGS 



Kennicutt et al.l 120031 ) . Overlaying the HI hole and iRhode et al.l (119991 ) aperture locations 



over these images, we compared the location of bright regions with known HI holes (Figures 
[29] — [3l]). Because the correlation between low SFRs and Ha, 24/xm, and UV is not well 
understood, we are only able to make qualitative conclusions. From all the images, it is 
evident that there is not a one to one correlation between any of these fluxes and the stellar 
populations inside the HI holes as shown in the HST/ACS based CMDs. This simple test 
demonstrates that deep photometry is necessary to understand the formation and histories 
of HI holes. 



5.4. Comparison with the SMC and LMC 



Perhaps the biggest challenge to the stellar feedback model of HI hole creation comes 
from observations of some of the closest galaxies. The proximity of the SMC and LMC 
allow for detailed observatio ns, including high resolution HI imaging, which probes physical 



15 pc (iKim et al.l Il999l ). Matching catalogs of young stellar objects 



scales as small as 

(e.g., OB associations, sup ergiants, WR st ars) and Ho ; imaging to HI hole catalo gs, three 



recent studies of th e LMC flKim et al.lfl999h. the SMC flHatzidimitriou et aLllioOsh . and the 
Magellanic Bridge (IMuUer et al.,.2003,) do not find strong correlations between the locations 
of young stellar objects and HI holes. The largest challenge to the stellar feedback creation 
hypothesis posed by these studies a re the presence of young HI holes (~ 2—10 Myr) with 
no associated yo ung stellar objects (IStanimirovig 120071 ). This is particularly problematic in 
the SMC, where iHatzidimitriou et al. I (l2005h identified 59 young holes (~ 5 — 10 Myr) that 



appear to be void of young stellar objects. IHatzidimitriou et al.l (120051 ) further find that 
the physical characteri stics of these particul a r hole s are no different than those with stellar 
components. Further, IHatzidimitriou et al.l (120051 ) conclude that the properties of the HI 
holes in the inner and outer region are indistinguishable and that from their analysis there 
is no evidence to suspect a different creation mechanism for inner and outer holes. 

Using the same type of analysis as in ^ 5.3, and the assumption that inner and outer holes 
are indistinguishable in their HI properties IHatzidimitriou et al.l (120051 ). we can compute the 
properties of the expected stellar populations if these young holes were formed due to stellar 
feedback. For the 59 empty HI holes in the SMC, we find that the mean size is ~ 128 
pc, the expansion velocity is 7.5 km s~^, and the inferred kinematic age is 8 Myr. We 
also compute Enoie ~ 56 xlO^°, assuming an HI volume density of 1.0 cm~^, and find that 
SFRxA = 7xl0~^ M0 yr~^. We simulate 100 CMDs assumi ng a distance modulus of 19, 
photometric completeness limits of my = 22 and m/ = 21.7 (IHarris fc Zaritskyll2004l ). and 
the parameters listed in §3.2. We chose to only consider a burst model of SF at 8 Myr, since 
the differences between a burst and constant model of SF over such a short time scale are 
minimal. We ran 100 simulations, to account for small number statistics, and found that the 
typical CMD contained ~ 40 stars (within the completeness limits) and the mean magnitude 
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of the brightest star in the CMD was uip^bbw ~ 14.80. We found that only 25 out of the 100 
simulated CMDs contained at least one star with vcipb^bw < 14, i.e., a bright OB type star. 
For TdFbbbW < 15, we found 72 CMDs contained at least one such star. This test signifies 
that for a episode of SF 8 Myr ago, only 25% of the time we would expect to find a star 
brighter than mp^bbw = 14, and 72% of the time for \iiFbbbW = 15. 



Applying these statistics to the 509 holes cataloged by iHatzidimitriou et al.l (120051 ) . we 
expect 128 HI holes with at l east one star brighte r than vcip^^^w = 14, or 367 holes for 
^F555W = 15. In comparison, IHatzidimitriou et al.l (l2005l ) find 330 HI holes with at least 
one star cluster, young star, or WR star. Initially, they also found 80 holes that appeared 
void of stars, however follow up observations revealed that 21 of the 80 holes appeared to 
have A spectral-type stars, leaving 59 holes without detected stellar companions, which is 
consistent with the 143 of holes we predict will appear void of MS stars at mp^^^^r = 15. The 
uncertainties in our calculation are at least 50%, due to uncertainties in Euoie, the stellar 
evolution models, and statistical error. Further, we are not able to test out results directly 
as we cannot reasonabl y compare our predicted ste llar populations to the various catalogs 
used for correlations by IHatzidimitriou et al.l (120051 ). However, we have demonstrated that 
even without the detection of a bright object inside a young HI hole, it is possible that the 
hole was created by stellar feedback. 



6. Conclusions 

In summary: 



The H ST/ ACS based CMDs of sta rs within HI holes from the catalogs of iPuche et al. 



(I1992I ) and iBagetakos et al.l (120091 ) revealed significant stellar populations, including 
young MS and BHeBs with ages less than ~ 50 Myr, in all holes. We found that ^ 1% 
of the stars in the composite stellar population of either holes catalog were luminous 
MS stars less than ~ 10 Myr in age. This indicates that young OB associations are 
not reliable tracers of the locations of HI holes. 

We searched for differences in the stellar populations inside HI holes and those in 
selected control fields, regions similar in size to holes, but that span a wide range 
in HI column density. We found that, on average, the control fields have twice as 
many luminous (i.e., young) MS stars as the HI holes, however this is consistent with 
statistical fluctuations from small number statistics tests. We measured the SFHs and 
cumulative SFHs of both samples (holes and control fields) and found no significant 
differences. However, because we are not sensitive to holes less than ~ 100 pc in size, 
it could be that smaller holes are already forming in control fields. Only further similar 
studies involving an even wider range of hole sizes (e.g., SMC or LMC) can address 
this question with the requisite level of precision. 
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Using the SFHs for stars in each hole as input into STARBURST99, we calculated the 
energy associated with SF over the past 200 Myr. From this we computed the energy 
input over the inferred kinematic age of each hole and calculated the feedback efficiency. 
Compared to theoretical values, 8 of the 23 P92 and 13 of the 19 THINGS holes have 
efficiencies greater 20%, with several exceeding 100%. Looking at the energy profiles 
over the last 200 Myr, we found that stellar feedback did produce enough energy, in 
principle, to have created all the HI holes, just not necessarily within the inferred 
kinematic age. This finding suggests that the inferred kinematic age (which is poorly 
constrained from observations) may not represent the true age of an individual HI hole. 
Indeed, the concept of an age for an HI hole created by multiple generations of SF is 
intrinsically ambiguous. 

Using the BHeBs from the full HST/ACS CMD of Ho II, we constructed two spatially 
resolved recent SFH movies of Ho II with look back times of ~ 100 Myr and ~ 200 Myr. 
Comparing these movies with the HI distribution, we found that low levels of SF have 
been common, and only in the past ~ 40 Myr has the SF rate elevated significantly. 
Since many of the holes are older than ~ 40 Myr, this suggests that holes are formed 
through a series of small SF events as opposed to a larger event at a single epoch. 

From th e HST/ACS photora etry, we constructed CMDs corresponding to the apertures 



used by iRhode et al.l (Il999l ). In all cases, we found mixed age stellar populations of 
hundreds or thousands of stars. Using Monte Carlo tests of synthetic CMDs and em- 
pirical comparisons from the HST/ACS photometry, we demonstrated that integrated 
light does not accurately trace the SFHs of stellar populations inside HI holes. We also 
show that young MS stars (< 75 Myr) are not highly clustered, which may imply that 
the single age stellar cluster met hod of HI formation is not the most likely model. This 



assumption can account for why iRhode et al.l (119991 ) did not detect a greater number 
of holes with stars inside. 



Using the Monte Carlo technique, we computed expected stellar populations for HI 
holes outside the optical body of Ho II and found that low levels of SF (~ 10~^ M0) 
could account for the presence of those holes. The remnant populations would be quite 
faint as only a handful of young stars would remain to be seen at the present time. We 
also examine the use of other SF tracers. Ha, 24^m, and UV, and find that at such 
low SFRs, they do not reliably correlate with the stellar populations of HI holes. 

We applied s imilar analysis to the 59 y oung (~ 8 Myr) HI holes that appear empty 



in the SMC flHatzidimitriou et all 120051 ). finding SFR^^a = 7x10"^ M© yr-\ Using 
the integrated stellar mass, we assume a burst of SFH at 8 Myr, and construct 100 
simulated CMDs. The resulting synthetic stellar populations show that only 25% of 
the CMDs have at least one MS stars brighter than mp^^^w = 14 and 72% have at 
least one MS star brighter than mF555w = 15. Applying this to the 509 holes in the 
sample, we expect 128 HI holes with at least one bright star in the 25% case, 367 holes 
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in the 72% case, and 143 ho les to be without any MS star brighter than m^sssvi/ = 15. 



Hatzidimitriou et aLl (120051 ) find 330 holes that have at least one bright object, and 
59 HI holes that do not have any objects. While this comparison is subject to large 
uncertainties, it demonstrates the plausibility of a stellar feedback origin for young HI 
holes without a bright stellar companion. 

The evidence presented by the HST/ACS observations of Ho II suggests that stellar 
feedback provides enough energy, in principle, to create the observed HI holes in Ho II. The 
traditional theory that single epoch episodes of SF, and singled aged clusters, are responsible 
for HI hole creation is unlikely, because of the presence of mixed age stellar popu lations inside 



all of the holes, including holes with younger inferred ages. The conclusions of iRhode et al. 



(119991 ) and those in this paper are not contradicting, rather they both seem to reinforce the 
idea that HI holes are not often created by single aged stellar populations. The summation 
of the evidence is that the likely mechanism for the creation of HI holes is stellar feedback 
from SNe due to multiple generations of SF. An important caveat is found in the similarity 
of the stellar populations in control fields and the HI holes. While this detracts from the 
potential casual relationship between stellar feedback and HI hole creation, this result could 
simply be explained by the fact that we are not sensitive to holes smaller than ~ 100 pc. 
Applying similar analysis to galaxies with a wider range in HI holes sizes (e.g., SMC or 
LMC) could reveal the extent of this effect and provide further insight into the potential link 
between stellar feedback and the formation and growth of HI holes. 
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Table 1. Photometric Properties of Puche et al. HI Holes 



Hole 


a 


5 


Diameter 


No. of Stars 


Integrated 


Number 


(J2000) 


(J2000) 


(arcsec) 


in CMD 


T^F555W 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


10 


8 18 42 


+70 41 21 


73 
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17.68 


12 


8 13 49 


+70 43 33 


22 


2065 


18.65 


14 


8 18 56 


+70 43 27 


27 


3907 


17.47 


16 


8 18 59 


+70 43 51 


13 


948 


19.00 


20 


8 19 01 


+70 42 44 


8 


383 


20.19 


01 




_i_7n /II ^ fi 
-r (yJ 41 ID 


lie: 
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Q01 c;q 


1 A 70 


22 


8 19 06 


+70 43 08 


32 


6443 


16.65 


23 


8 19 06 


+70 43 30 


39 


9477 


16.22 


24 


8 19 08 


+70 44 03 


10 
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19.37 


28 


8 19 13 


+70 43 31 


19 
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17.93 


29 


8 19 14 


+70 42 21 


27 


4517 


17.58 


30 


8 19 14 


+70 44 20 


69 


20614 


15.77 


31 


8 19 15 


+70 43 49 


18 


1625 


18.37 


32 


8 19 16 


+70 41 49 


37 


3782 


17.83 


33 


8 19 17 


+70 44 56 


25 


2140 


18.34 


36 


8 19 19 


+70 43 14 


37 


5945 


16.98 


37 


8 19 20 


+70 45 17 


20 


1011 


19.12 


39 


8 19 23 


+70 44 43 


23 


1971 


18.24 


42 


8 19 25 


+70 45 12 


22 


1383 


18.75 


43 


8 19 27 


+70 42 06 


16 


549 


18.75 


44 


8 19 27 


+70 43 31 


48 


7566 


16.82 


45 


8 19 27 


+70 45 19 


32 


2344 


18.08 


48 


8 19 36 


+70 44 54 


50 


4695 


16.93 



Note. — Numbe rs, coordinates, and diameters for each hole taken from 
Puche et al. (|l992h . Hole coordinates have been converted to the J2000 
epoch for ease of comparison. 
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Table 2. Photometric Properties of THINGS HI Holes 



Hole 


a 


S 


Diameter 


No. of Stars 


Integrated 


Number 


(J2000) 


(J2000) 


(arcsec) 


in CMD 


^F555W 


(1) 


(2) 
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(5) 
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+70 41 36 


44 


1733 


18.81 


6 


8 18 49 


+70 43 40 


29 


3781 


17.98 
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19 
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18.62 


31 
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18.05 


34 


8 19 35 


+70 45 05 


40 


3673 


17.33 


39 


8 19 44 


+70 43 55 


29 
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19.31 



Note. — HI holes from THINGS (|Walter et al.ll2008l ) that are within 
the HST/ACS field of view. Numbers, co ordinates, and diameters for each 
HI hole are from iBagetakos et al. ( 20091) 



Table 3. Derived Properties of Puche et al. HI Holes 



Hole Diameter DV Kinematic Age ^Hole {P92) T^Hole (new) T^SFKA Putative 
Number (pc) km s^^ (Myr) {10^° erg) (lO'^" erg) (lO^^" erg) Efficiency^A 

(1) (2) (3) (4) (5) (6) (7) (8) 



10 


1171 


8.8 


65 


500 


364 


2400 


15% 


12 


349 


7.7 


22 


18.0 


16 


130 


12% 


14 


421 


3.9 


53 


15.0 


19 


2700 


0.7% 


16 


199 


1.0 


97 


0.30 


0.5 


1000 


<0.1% 


20 


125 


6.5? 


9.4 


0.4 


2.0 


10 


20% 


21 


1820 


6.6 


135 


2800 


7700 


60000 


13% 


22 


508 


6.5 


38 


48.0 


230 


1600 


15% 


23 


619 


4.5 


67 


60 


211 


3300 


6.4% 


24 


162 


7.5? 


11 


1.66 


5.0 


20 


27% 


28 


302 


7.2? 


21 


11.0 


83 


220 


37% 


29 


469 


3.9 


59 


21.4 


190 


1300 


15% 


30 


1083 


8.1 


65 


737 


2200 


5700 


39% 


31 


289 


2.5? 


56 


2.20 


13 


900 


1.4% 


32 


583 


4.8 


59 


48.4 


420 


2600 


16% 


33 


388 


3.5 


54 


10.2 


17 


780 


2.2% 


36 


579 


8.3 


34 


119 


990 


1800 


55% 


37 


312 


11.9 


13 


28.8 


34 


45 


74% 


39 


359 


3.5 


50 


8.00 


16 


1500 


1.1% 


42 


355 


10.1 


17 


34.0 


39 


75 


52% 


43 


255 


4.8 


26 


3.73 


25 


90 


28% 


44 


752 


7.9 


47 


217 


1200 


7300 


17% 


45 


596 


6.6 


37 


53.5 


51 


930 


6% 


48 


786 


14.1 


27 


649 


560 


450 


120% 



Note. — (1) Hole numbers from [Puche et al.l lll992 
Mpc; (3) Expansion velocities from iPuche et al.l 1 1992 



(2) Hole diameters calculated assuming a distance of 3.4 
; (4) Kinemati c age assumes the diameter from column 2 



and expansion velocity from column 3; (5) Energy from [Puche et al.l l|l992f ) was adjusted for a distance of 3.39 
Mpc; (6) Sffoie computed using Equation 1 and the diameters and expansion velocities in this table, along with 
the values of no computed in §4.1; (7) Total energy over the inferred kinematic age of the HI hole computed using 
STARBURST99; (8) The ratio of EhoU to Eki„^^atic- 
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Table 4. Derived Properties of THINGS HI Holes 



Hole 


Diameter 


DV 


Kinematic Age 


EhoIc (P92) 


^Hole (new) 


^SFKA 


Putative 


Type 


Number 


(pc) 


km s^^ 


(Myr) 


(IqSO erg) 


(IqSO erg) 


(IqSO erg) 


Efficiency^ A 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


4 


731 


7 


51 


100 


170 


500 


34% 


1 


6 


486 


7 


34 


50 


180 


600 


31% 


1 


8 


345 


10 


17 


32 


140 


500 


28% 


3 


9 


279 


10 


14 


13 


54 


15 


360% 


3 


12 


403 


11 


18 


63 


450 


330 


140% 


3 


14 


318 


15 


10 


40 


310 


20 


1600% 


2 


16 


761 


7 


53 


200 


2100 


4600 


45% 


1 


17 


2110 


7 


147 


3200 


18000 


100000 


18% 


1 


19 


1107 


7 


77 


500 


3200 


9500 


34% 


1 


20 


329 


12 


13 


25 


56 


100 


56% 


3 


21 


378 


9 


21 


20 


88 


300 


29% 


3 


23 


738 


16 


23 


630 


2900 


2400 


120% 


2 


26 


483 


7 


34 


50 


133 


2200 


6% 


1 


27 


312 


7 


22 


20 


25 


190 


13% 


1 


29 


479 


7 


33 


32 


43 


900 


5% 


1 


30 


318 


13 


12 


32 


40 


150 


26% 


2 


31 


491 


7 


34 


50 


83 


1000 


8% 


1 


34 


657 


7 


46 


80 


79 


730 


11% 


1 


39 


474 


18 


13 


160 


66 


50 


130% 


2 



Note. — (1) — (5) Hole number, size, expansion, velocity, and Effole from the catalog of iBagetakos et al.l ||2009| ): (6) 
Total energy over the inferred kinematic age of the HI hole computed using STARBURST99; hole 6 had no appreciable SF 
within its kinematic age; (6) -Enoie computed using Equation 1 and the diameters and expansion velocities in this table, 
along with the values of no computed in §4.1; (7) Total energy over the inferred kinematic age of the HI hole computed 
using STARBURST99; (8) Feeback efficiency: The ra tio of Enoie to Esfka ; (8) The ratio of SFRka to the SFR averaged 
over the past 500 Myr. (9) Hole type as classified bv IBagetakos et al.l | |2009|) . where type 1 holes arc blown out and type 2 
or 3 holes are still appreciably expanding. 
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Table 5. Photometric Properties of Control Fields 



Hole 


a 


S 


Aperture Diameter 


No. of Stars 


Integrated 


Number 


(J2000) 


(J2000) 


(arcsec) 


in CMD 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


cl 


8 19 46 


+70 44 29 


30 


839 


19.38 


c2 


8 19 39 


+70 43 55 


30 


1309 


19.20 


c3 


8 19 28 


+70 44 09 


30 


4089 


17.64 


c4 


8 19 35 


+70 43 08 


30 


1549 


18.78 


c5 


8 19 32 


+70 42 40 


30 


1863 


18.08 


c6 


8 19 21 


+70 42 28 


30 


3598 


17.96 


c7 


8 19 12 


+70 42 55 


30 


5183 


17.40 


c8 


8 18 51 


+70 42 05 


30 


1988 


18.72 


c9 


8 18 45 


+70 42 39 


30 


2147 


17.74 



Note. — Control fields shown in cyan in Figure [21 Corresponding CMDs in Figure 
[Hand SFHs in Figure [HI 



Table 6. Photometric Properties of HI Holes used by lRhode et al.l (119991 ) 



Hole a b Aperture Diameter No. of Stars Hole Integrated 

Number (J2000) (J2000) (arcsec) in CMD Type mpsssw 

(1) (2) (3) (4) (5) (6) (7) 



10 


8 18 


40 


+70 


41 


39 


30.00 


893 


a 


19.75 


12 


8 18 


50 


+70 


43 


41 


20.00 


1989 


b 


18.63 


14 


8 18 


55 


+70 


43 


34 


20.00 


2810 


c 


18.12 


16 


8 18 


58 


+70 


43 


01 


20.00 


2740 


d 


17.12 


20 


8 18 


59 


+70 


42 


52 


20.00 


2639 


c 


17.44 


21 


8 19 


02 


+70 


41 


21 


30.00 


4077 


d 


17.18 


22 


8 19 


04 


+70 


43 


16 


20.00 


3064 


d 


17.24 


23 


8 19 


04 


+70 


43 


40 


20.00 


3051 


c 


17.85 


24 


8 18 


07 


+70 


44 


14 


20.00 


2825 


c 


18.08 


28 


8 19 


12 


+70 


43 


52 


20.00 


2526 


c 


18.01 


30 


8 19 


13 


+70 


44 


32 


20.00 


2062 


c 


18.59 


31 


8 19 


14 


+70 


44 


04 


20.00 


2374 


c 


18.09 


32 


8 19 


12 


+70 


41 


54 


24.00 


2615 


c 


18.31 


33 


8 19 


16 


+70 


45 


10 


20.00 


1628 


c 


18.81 


36 


8 19 


19 


+70 


43 


21 


24.00 


3324 


d 


17.65 


37 


8 19 


19 


+70 


45 


32 


20.00 


1205 


b 


19.04 


39 


8 19 


22 


+70 


44 


51 


20.00 


1893 


c 


18.29 


42 


8 19 


21 


+70 


45 


21 


20.00 


1336 


b 


18.76 


43 


8 19 


26 


+70 


42 


16 


16.00 


805 


8 


19.40 


44 


8 19 


24 


+70 


43 


33 


30.00 


4229 


d 


17.60 


45 


8 19 


25 


+70 


45 


37 


20.00 


1023 


c 


19.04 


48 


8 19 


35 


+70 


45 


06 


30.00 


1892 


d 


18.14 



Note. — ( 1) — (4) Numbers, coordinates, and diameters for each apertur e take n from 
Rhode et al. 1 (|l999l) . which are based on the HI hole catalog of IPuche et all (Il992l) . Co- 
ordinates have bee n converted to J2000 for ease of comparison; (6) Classifications from 



Rhode et al.l (|l999[ ) 



— empty hole, b — faint foreground star, c — galaxian background, 
d — possible star cluster, e — possible photoionization region; (7) integrated magnitudes 
calculated from the ACS CMDs (Figures [25] & [ 
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Fig. 1. — (a) HI surface d ensity image of Ho II processed by the THINGS program 
(THINGS; IWalter et all l2008h with a contour of IQ-^^ cm'^ overlaid; ( b) HI holes cat 
alogu ed by iPuche et al.l (119921): fc) HI holes fro m the THINGS catalog (IBagetakos et al. 
20091 ): (d) apertures used by iRhode et al.l ( 1l999l ) with the ACS footprint overlaid in black. 
Note the 'shredded' appearance of the HI distribution as well as the numerous separate holes 
cataloged in regions of low HI column density. 
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Fig. 2. 

10-21 



(a) THINGS processed HI image of Ho II flWalter et al.l 120081) wi t h a c ontour of 
cm '"^ overlaid and (b) ACS drizzled image of Ho II with iPuche et al.l (Il992[) holes in 
blue and control fields in cyan; (c) and (d) THING hole s catal og (IBagetakos et al.ll2009l ) in 
green and control fields in cyan; (e) and (f ) iRhode et al.l (119991 ) apertures in red and control 
fields in cyan. 
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Fig. 3.— HST/ACS CMDs of the stars inside th e HI holes cataloged bv lPuche et all (119921 ) . 
corrected for foreground reddening, Ay = 0.11 (ISchlegel et al.lll998l ). The schematic in the 
upper right hand corner shows the ages of each type of star, MS (green), BHeB (green), and 
RHeB (red), with the logarithm of the ages (turn off age for the MS) shown as a function of 
magnitude and color on the CMD. Although sparse in some fields, note the presence of young 
MS and BHeB stars in all the CMDs. The fact the the BHeBs span a range in magnitudes 
indicate that multiple episodes of recent SF must have taken place, as different age BHeBs 
do not overlap on the CMD. In contrast, BHeBs from a single age cluster would have an 
overdensity at only one magnitude. 
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Fig. 4.— HST/ACS CMDs of the stars inside th e HI holes cataloged bv lPuche et all (11992h . 
corrected for foreground reddening, Ay = 0.11 (jSchlegel et al.l 119981 ). The schematic in the 
upper right hand corner shows the ages of each type of star, MS (green), BHeB (green), and 
RHeB (red), with the logarithm of the ages (turn off age for the MS) shown as a function of 
magnitude and color on the CMD. Although sparse in some fields, note the presence of young 
MS and BHeB stars in all the CMDs. The fact the the BHeBs span a range in magnitudes 
indicate that multiple episodes of recent SF must have taken place, as different age BHeBs 
do not overlap on the CMD. In contrast, BHeBs from a single age cluster would have an 
overdensity at only one magnitude. 
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Fig. 5. — HST/ACS CMDs of the stars inside the HI holes cataloged by iBagetakos et al. 



( 2009 ). corrected for foreground reddening, Ay = 0.11 ( Schlegel et al. 19981 ). The schematic 
in the upper right hand corner shows the ages of each type of star, MS (green), BHeB 
(green), and RHeB (red), with the logarithm of the ages (turn off age for the MS) shown 
as a function of magnitude and color on the CMD. Although sparse in some fields, note 
the presence of young MS and BHeB stars in all the CMDs. The fact the the BHeBs span 
a range in magnitudes indicate that multiple episodes of recent SF must have taken place, 
as different age BHeBs do not overlap on the CMD. In contrast, BHeBs from a single age 
cluster would have an overdensity at only one magnitude. 
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Fig. 6.— HST/ACS CMDs of the nine control fields in Ho II (T able [3] and Figure [2]), 
corrected for foreground reddening, Ay = 0.11 (jSchlegel et al.lll998l ). Although sparse in 
some fields, note the presence of young MS and BHeB stars in all the CMDs. The fact the 
the BHeBs span a range in magnitudes indicate that multiple episodes of recent SF must 
have taken place, as different age BHeBs do not overlap on the CMD. In contrast, BHeBs 
from a single age cluster would have an overdensity at only one magnitude. Comparing the 
CMDs of the control fields to those associated with the HI holes (Figure [3] — [5]), there is 
not a clear difference in the composition of the two stellar populations despite having vastly 
different HI column densities. 
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Fig. 7.— HST/A CS CMD of THINGS hole 17 with the Z=0.002 isochrones from 



Marigo et al.l (120081 ) overlaid. The observed CMD has been corrected for foreground red- 



dening. Notice that although the observed MS and isochrones show excellent agreement, the 
BHeB models become increasingly redder with age, when compared with the data. This is 
not an effect of rnetalli city, as we matched the isochrones to the measured nebular abundance 
(IMiller fc Hodgelll996l ). Instead, this is an example of the color mismatch for BHeB models 
at metallicites greater than ~ 10% Z©. 
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Fig. 8. — Comparison of the observed (panel (a)) and model CMDs (panel (b)) for THINGS 
hole 23 as measured using the code of iDolphinI ([20021). The bottom panels show the residual 
CMD (panel (c)) and the significance of the residuals (panel (d)). In the bottom panels 
black points indicate many more observed stars than synthetic star, while white points have 
more synthetic than observed stars. From the residual significance diagram we do not see 
any particularly poorly matched areas, indicating the synthetic CMD matches the data quite 
well. 
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Fig. 9. — The recent SFHs of the stars inside the HI holes cataloged by lPuche et al.l (1l992l ) 
over the last 200 Myr. The time resolution is ~ 10 Myr until 50 Myr ago, and then ~ 25 
Myr there after. The red dotted line indicates the inferred kinematic age of each hole (Table 
[3]). Note that the SFH of hole 21 has been scaled down by a factor of 10. In most cases, the 
majority of SF has taken place more recently than the inferred kinematic age, however, this 
is not always the case. 
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Fig. 10. — The recent SFHs of the stars inside the HI holes cataloged by lPuche et al.l ( 1l992l ) 
over the last 200 Myr. The time resolution is ~ 10 Myr until 50 Myr ago, and then ~ 25 Myr 
there after. The red dotted line indicates the inferred kinematic age of each hole (Table [3]). 
In most cases, the majority of SF has taken place more recently than the inferred kinematic 
age, however, this is not always the case. 
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Fig. 11. — The recent SFHs of the stars inside the HI holes cataloged by iBagetakos et al. 
(120091) over the last 200 Myr. The time resolution is ~ 10 Myr until 50 Myr ago, and then 
~ 25 Myr there after. The red dotted line indicates the kinematic age of each hole (Table 
H]). Note that the SFH of hole 17 has been scaled down by a factor of 10. In most cases, the 
majority of SF has taken place more recently than the inferred kinematic age, however, this 
is not always the case. 
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Fig. 12. — The recent SFHs of each control field in Ho II over the last 200 Myr. The time 
resolution is ~10 Myr ber bin until 50 Myr and ~ 25 Myr there after. 
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Fig. 13. — Compar ison of the cumulative SF of stars within t he holes from the Ho I I for the 
Puche et all fll992f ) catalog (blue) and the THINGS catalog fiBagetakos et alJl2009f ) (green) 



with the control fields cyan. The left panels show the cumulative SF over the the past 14.1 
Gyr, while the panels on the right show the cumulative SF over the past 200 Myr. Although 
any given cumulative SFH is unique, a comparison of the HI holes and control fields suggest 
that the SFHs of the stars inside HI holes and the control fields are indistinguishable. For 
reference, the globally averaged SFR for the galaxy is 3. 6xl0~^ Mq yr~^, and the full SFH 
of the entire galaxy was measured by lWeisz et al.l (l2008l ). 
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Fig. 14. — Top panel: Surface brightness pro file of Ho II as measured from 3 . 6 Mm imaging 
from the Local Volume Legacy program (LVL: lLee et al.ll2008l : iDale et al.ll2009l ) as a function 
of radius; Second Panel: The mass surfac e density derived from the surface brightness profile 
following the method of lOh et al.l (120081 ). The red points are the data points and the green 
line is the best fit analytic function; Third panel: The gas scale hei ght of Ho II as a function 
of rad ius assuming an isothermal disk derived using the method of ( Kellman 1970l : Kim et al. 



19991 ). Data points are in red and the best fit analytic f unction is the gree n line. The grey 



dotted line is the constant value gas scale height used by ( Puche et al. 19921 ): Bottom Panel: 
The pre sent day HI yolurn e density, function of radius (green line) with the values 

used by iPuche et al.l (119921 ) in magenta. Note the difference in the valu es of n between the 
best fit model obtained via surface photometry and the values used by iPuche et al.l (119921 ) 
in the inner ~ 2 kpc. Because EhoIs oc /Zq^^ (where uq is the average volume density prior 
to hole formation), an increase in the volume density leads to an increase in the minimum 
energy necessary to create an HI hole. In this case no measures the average value of the HI 
volume density, however, SF is more likely to occur at HI peaks, thus Enoie is likely a lower 
bound (See §4.1 for more details). 
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Fig. 15. — The cumulative energi es of the stars inside the HI holes of iPuche et al.l (Il992l ) 



computed using STARBURST99 flLeitherer et al. 1999 ) with the measured SFHs as input. 
The energies are integrated from the present (t = 0) over the past 200 Myr. The vertical grey 
dashed line is the kinematic age of each hole. The horizontal dotted lines are Ehqi^ assuming 
100% (red), 10% (green), and 1% (blue) efficiencies. The intersection of the grey line with 
the energy profile indicates the putative feedback efficiency over the inferred kinematic age. 
In several cases, this falls outside the expected range, and may even be greater than 100%. 
This suggests that inferred kinematic age may not represent the true age of the hole. If 
instead, we use the intersection 10% efficiency line (green) and the energy profile as a guide, 
we can get a rough sense of how old the hole may be, if it formed from multiple generations 
of SNe. 
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Fig. 16. — The cumulative energi es of the stars inside the HI holes of iPuche et al.l (Il992l ) 
computed using STARBURST99 fiLeitherer et al. 1999 ) with the measured SFHs as input. 
The energies are integrated from the present (t = 0) over the past 200 Myr. The vertical grey 
dashed line is the kinematic age of each hole. The horizontal dotted lines are Ehqi^ assuming 
100% (red), 10% (green), and 1% (blue) efficiencies. The intersection of the grey line with 
the energy profile indicates the putative feedback efficiency over the inferred kinematic age. 
In several cases, this falls outside the expected range, and may even be greater than 100%. 
This suggests that inferred kinematic age may not represent the true age of the hole. If 
instead, we use the intersection 10% efficiency line (green) and the energy profile as a guide, 
we can get a rough sense of how old the hole may be, if it formed from multiple generations 
of SNe. 



- 51 - 




Fig. 17. — The cumulative energi es of the stars inside the HI holes of iPuche et al.l (Il992l ) 



computed using STARBURST99 fiLeitherer et al. 1999) with the measured SFHs as input. 
The energies are integrated from the present (t = 0) over the past 200 Myr. The vertical grey 
dashed line is the kinematic age of each hole. The horizontal dotted lines are Ehqi^ assuming 
100% (red), 10% (green), and 1% (blue) efficiencies. The intersection of the grey line with 
the energy profile indicates the putative feedback efficiency over the inferred kinematic age. 
In several cases, this falls outside the expected range, and may even be greater than 100%. 
This suggests that inferred kinematic age may not represent the true age of the hole. If 
instead, we use the intersection 10% efficiency line (green) and the energy profile as a guide, 
we can get a rough sense of how old the hole may be, if it formed from multiple generations 
of SNe. 
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Fig. 18. — The cumulative energie s of the stars inside th e HI holes of iBagetakos et al.l (120091 ) 
computed using STARBURST99 flLeitherer et al. 1999 ) with the measured SFHs as input. 
The energies are integrated from the present (t = 0) over the past 200 Myr. The vertical grey 
dashed line is the kinematic age of each hole. The horizontal dotted lines are Ehqi^ assuming 
100% (red), 10% (green), and 1% (blue) efficiencies. The intersection of the grey line with 
the energy profile indicates the putative feedback efficiency over the inferred kinematic age. 
In several cases, this falls outside the expected range, and may even be greater than 100%. 
This suggests that inferred kinematic age may not represent the true age of the hole. If 
instead, we use the intersection 10% efficiency line (green) and the energy profile as a guide, 
we can get a rough sense of how old the hole may be, if it formed from multiple generations 
of SNe. 
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Fig. 19. — The cumulative energie s of the stars inside th e HI holes of iBagetakos et al.l (120091 ) 
computed using STARBURST99 flLeitherer et al. 1999 ) with the measured SFHs as input. 
The energies are integrated from the present (t = 0) over the past 200 Myr. The vertical grey 
dashed line is the kinematic age of each hole. The horizontal dotted lines are Ehqi^ assuming 
100% (red), 10% (green), and 1% (blue) efficiencies. The intersection of the grey line with 
the energy profile indicates the putative feedback efficiency over the inferred kinematic age. 
In several cases, this falls outside the expected range, and may even be greater than 100%. 
This suggests that inferred kinematic age may not represent the true age of the hole. If 
instead, we use the intersection 10% efficiency line (green) and the energy profile as a guide, 
we can get a rough sense of how old the hole may be, if it formed from multiple generations 
of SNe. 
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Fig. 20. — The cumulative energie s of the stars inside th e HI holes of iBagetakos et al.l (120091 ) 
computed using STARBURST99 flLeitherer et al. 1999 ) with the measured SFHs as input. 
The energies are integrated from the present (t = 0) over the past 200 Myr. The vertical grey 
dashed line is the kinematic age of each hole. The horizontal dotted lines are Ehqi^ assuming 
100% (red), 10% (green), and 1% (blue) efficiencies. The intersection of the grey line with 
the energy profile indicates the putative feedback efficiency over the inferred kinematic age. 
In several cases, this falls outside the expected range, and may even be greater than 100%. 
This suggests that inferred kinematic age may not represent the true age of the hole. If 
instead, we use the intersection 10% efficiency line (green) and the energy profile as a guide, 
we can get a rough sense of how old the hole may be, if it formed from multiple generations 
of SNe. 
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Fig. 21. — Putative stellar feedback efficiencies {Esfka/ Euoie] Tables |3] and H]) plotted versus 
the galactrocentric radius (Rgc); h ole radius i^Hnip), H I expansion velocity iy Expansion) , and 
kinematic age of each hole for the iPuche et al.l and the THINGS HI hole catalogs. 

For the THINGS holes, green stars are type 1 holes (blown out), the red stars are type 2 
or 3 (not blown out), and the arrow indicate holes that have efficiencies greater than 100% 
and have been moved to fit the plot scale. The vertical grey dashed lines i i idicate the range 



of efficiencies predicted by various mode ls, 1% — 20% (ITheis et al.lll992l : ICole et al.lll994 



Padoan et al.l 1 199 71 : iThornton et al.lll998f ). Note that the type 2 and 3 holes in the THINGS 



hole sample do not fall in to the expected range of putative feedback efficiency. These holes 
also tend to be young, expanding quickly, and are generally smaller in size. Arrows indicate 
lower limits on the efficiency for both hole catalogs. 
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Fig. 22. — The ra tio of look back age to infer red kinematic age for the HI holes from the 



Puche et al.l ( 1l992l ) and iBagetakos et al.l (120091) as a function of the inferred kinematic age. 



The look back age is defined as the time at which the line of 10% efficiency intersects the 
energy profiles in Figures [15] — [201 The spread of values greater than unity (grey dashed line) 
implies that the inferred kinematic age often underestimates the age of HI holes, suggesting 
that the concept of an age for an HI hole created by multiple generations of SF is intrinsically 
ambiguous. 
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Fig. 23. — The spatially resolved recent SFH of Ho II over the last 90 Myr with 10 Myr time 
resolution as seen on the sky in units of SFR/area. The HI 10~^^ cm~^ contour is overlaid in 
white. Note that majority of the galaxy is dominated by constant low level SF, while only 
a few regions have experienced elevated episodes of SF. As expected, the most recent frame 
shows the peaks of SF coincident with the peaks of the HI distribution. 
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Fig. 24. — The spatially resolved recent SFH of Ho II over the last 180 Myr with 20 Myr 
time resolution as seen on the sky in units of SFR/area. The HI 10~^^ cm~^ contour is 
overlaid in white. Note that majority of the galaxy is dominated by constant low level SF, 
while only a few regions have experienced elevated episodes of SF. As expected, the most 
recent frame shows the peaks of SF coincident with the peaks of the HI distribution. 
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Fig. 25. — HST/ACS CMDs of the photomet ric apertures used b y iRhode et al.l (1l999l ). 
corrected for foreground reddening, Ay = 0.11 (ISchlegel et al.lll998l ). The schematic in the 
upper right hand corner shows the ages of each type of star, MS (green), BHeB (green), and 
RHeB (red), with the logarithm of the ages shown as a function of magnitude. Ahhough 
sparse in some fields, note the presence of young MS and BHeB stars in all the CMDs. The 
fact the the BHeBs span a range in magnitudes indicate that multiple episodes of recent 
SF must have taken place, as different age BHeBs do not overlap on the CMD. In contrast, 
BHeBs from a single age cluster would have an overdensity at only one magnitude. 
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Fig. 26. — HST/ACS CMDs of the photomet ric apertures used b y iRhode et al.l (1l999l ). 
corrected for foreground reddening, Ay = 0.11 (ISchlegel et al.lll998l ). The schematic in the 
upper right hand corner shows the ages of each type of star, MS (green), BHeB (green), and 
RHeB (red), with the logarithm of the ages shown as a function of magnitude. Ahhough 
sparse in some fields, note the presence of young MS and BHeB stars in all the CMDs. The 
fact the the BHeBs span a range in magnitudes indicate that multiple episodes of recent 
SF must have taken place, as different age BHeBs do not overlap on the CMD. In contrast, 
BHeBs from a single age cluster would have an overdensity at only one magnitude. 
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Fig. 27. — A histogram of the integrated ACS mp^^^w magnitudes from the HST/ACS CMDs 
corresponding to the photometric a pertures used by iRhode et aL The red shaded 



region represents the apertures that IRhode et al.l (119991 ) classify as containing putative star 
clusters (Table [6]). Note that they regions containing putative stellar clusters are among the 
brightest in the sample. 
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Fig. 28. — The spatial dis tribution of young MS stars (red < 10 Myr, blue < 75 Myr) plotted 
over (a) HI holes from t he lPuche et al.l (1l992l ) catalog; (b) HI holes from th e THINGS catalog 
( Bagetakos et al. 20091 ): (c) photometric apertures of [Rhode et al. ( 1999 ): Panel (d) is the 
CMD of the entire galaxy of Ho II with the young MS stars color-coded to match the 
spatial distributions with MS stars < 10 Myr in age in red and < 75 Myr in blue. Note the 
tendency of the red points to appear highly clustered, while the blue points are more widely 
distributed. Further, note that most of the HI holes do not contain highly clustered regions 
of young stars, rather these are found near the edges of the holes or near HI peaks. 
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Fig. 29. — HI holes of iPuche et al.l f|l992h fblue) and control fields (cy an) overlaid on Ho II 
images of (a) GALEX NUV flOil de Paz et al.l [2007): ( b) LVL Ha (jlee et aLl Eooi): (c) 
SINGS 24^m iKennicutt et all (l2003h : (d) THINGS HI flWalter et al.ll2008h . Despite that 
fact that all the HI holes contain hundreds or thousands of stars, from these images alone, it 
would seem that the HI holes are void of stars. At low SFRs (< 10~^) Mq yr~^, it seems that 
these tracers of recent SF no longer reliably correlate with the underlying stellar populations 
(see §5.3 for further discussion). 
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Fig. 30. — HI holes of iBagetakos et al.l (|2009h f green) and control fields ( cyan) overlaid o n 
Ho II images of (a) GALEX NUV dOil de Paz et al.l boOT i): (b) LVL Ha (lLp_eLaD l2Q08|) : 
(c) SINGS 24//m lKennicutt et all fl2003h : (d) THINGS HI flWalter et al.ll2008h . Despite that 
fact that all the HI holes contain hundreds or thousands of stars, from these images alone, it 
would seem that the HI holes are void of stars. At low SFRs (< 10~^) Mq yr~^, it seems that 
these tracers of recent SF no longer reliably correlate with the underlying stellar populations 
(see §5.3 for further discussion). 
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Fig. 31. — Photometric apertures of lRhode et al.l f|l999l) fred) and control fields (c yan') over- 
laid o n Ho II images of (a) GALEX NUV dOil de Paz et a llbnOTh : (h) LVL Ha Eee_eLaL 
2008h : (c) SINGS 24^m iKennicutt et all fl2003h : (d) THINGS HI flwalter et al.ll2008h . De- 
spite that fact that all the HI holes contain hundreds or thousands of stars, from these images 
alone, it would seem that the HI holes are void of stars. At low SFRs (< 10~^) Mq yr~^, it 
seems that these tracers of recent SF no longer reliably correlate with the underlying stellar 
populations (see §5.3 for further discussion). 
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